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NOMENCLATURE 


a  function  defined  by  eq  11a 
Ae  function  defined  by  eq  27f 

A)  constant  where  i  =  0, 1, ...  3 

Ai  function  defined  by  eq  25d  and  25e,  where  t  =  0, 1 

A  constant 

b  function  defined  by  eq  11b 
bi  constant  where  t  =  1^ 

B)  ith  constituent  of  the  mixture.  Subscripts  t  =  1, 2,  and  3  are  used  to  denote  uidrozen  water,  ice 
and  soil  minerals,  respectively 

Ci  heat  capacity  of  the  ith  constituent 
d  uiut  of  time  (day) 

d)  density  of  the  tih  constituent 

e  void  ratio 

fi  mass  flux  of  the  tth  constituent  relative  to  that  of  soil  minerals  where  t  s  1, 2 

/lo  mass  flux  of  water  in  the  unfrozen  part  of  the  soil 
I  function  defined  by  eq  16b 

k  thermal  conductivity  of  a  frozen  fiinge  defined  by  eq  9a 

kfj  thermal  conductivity  of  die  unfrozen  part  of  the  soil 

ki  thermal  conductivity  of  an  ke  layer 

kfn  limiting  value  of  k  defined  by  eq  9c 

Kf,  hydraulic  conductivity  in  the  unfrozen  part  of  die  soil 

fCj  empiriced  function  defined  by  eq  4a  where  i  =  1, 2 

Kii  limiting  value  of  as  x  approaches  while  x  is  in  Rj,  i  =  1,2 

Kjg  limiting  value  of  as  x  approaches  fto  while  x  is  in  R],  i  =  1,2 

L  latent  heat  of  fusion  of  water,  334  J  g-^ 

m  location  of  the  free  end  of  the  column 

Mi  name  of  a  model  where  t  =  1, 2, 3 

n  boundary  in  Rg 

Wj  boundary  with  i  =  0, 1  where  no  denotes  the  boundary  where  T  =  0®C  and  ttj  the  interface 
between  an  ice  layer  and  a  frozen  fringe 
nio  boundary  between  R^g  and  R^ 

Po  gravity  term,  0.098  [kPa/cm] 

pressure  of  the  tth  constituent  where  i  =  1, 2 
Pio  value  of  Pi  at  Hq 
Pi*  value  of  Pi  at  n 
^21  value  of  P2  at  ni 

iv 


r  rate  of  heave 
Rq  unfrozen  part  of  the  soil 
R]  frozen  fringe 

Rio  part  of  the  soil  bounded  by  Hq  and  rtio 
Rii  essential  frozen  fringe 
Rj  ice  layer 

R„  region  in  the  diagram  of  temperature  gradients  where  an  ice  layer  melts 

R,  region  in  the  diagram  of  temperature  gradients  where  the  steady  growth  of  an  ice  layer  occurs 

R»  boundary  between  R,  and  Ru 

R^*  boundary  between  Ro,  and  R, 

Ru  region  in  the  diagram  of  temperature  gradients  where  the  steady  growth  of  an  ice  layer  does 
not  occur 

S  property  of  a  given  soil 

SPq  segregation  potential  defined  by  eq  2a 

t  time 

T  temperature 

Ti  temperature  at  ni 

T,o  defined  by  eq  41 

calculated  value  of  Tj  from  the  measured  temperature  profile  in  R^ 

T 1  empirically  determined  value  of  T* 

T**  temperature  at  ni  when  eq  1  holds  true 
Tq  constant 

(T),  average  temperature  gradient  in  R] 

AT  defined  by  eq  3c 
X  spatial  coordinate 
y  variable  defined  by  eq  17c 
y  variable  defined  by  eq  38b 

y,  variable  defined  by  eq  3a 

a(0  trajectory  [ai(f),  ao(f)]  in  the  diagram  of  temperature  gradients 
Oq  absolute  value  of  the  temperature  gradient  at  ng 

tti  absolute  value  of  the  limiting  temperature  gradient  as  x  approaches  rt}  while  x  is  in  R2,  defined 

by  eq6 

ctf  absolute  value  of  the  temperature  gradient  near  »]  in  R2 
(Xu  absolute  value  of  the  temperature  gradient  near  ?io  in  Rg 

Po  defined  by  eq  11c 

y  constant,  1.12  MPa  “C"* 

Yi  defined  by  eq  24b 

5  thickness  of  a  frozen  fringe 

defined  by  eq  27b 

thickness  of  an  essential  frozen  fringe  defined  by  eq  27a 
6q  defined  by  eq  13c 


n  defined  by  eq  9b 
p  composition  of  the  soil 
Pi  bulk  density  of  the  tth  ccmstituent 
o  effective  pressure  defined  by  eq  1  and  16c 
^  empirical  function  of  T  defined  by  eq  14a 

^  valueof^atTs  T] 

^1  empirical  function  of  T  defined  by  eq  14b 

^11  value  of  at  T  =  Tj 
<|y,  variable  defined  by  eq  21b 
<Pi  variable  defined  by  eq  25f 

(Oj  dimeiuionless  quantity  defined  by  eq.lSa  through  18d  where  i  -  0,  3 

i  subscript  denotes  the  tth  coi\stituent  of  the  mixture  consisting  of  unfrozen  water  (i » 1),  ice  (i 
=  2)  and  soil  minerals  (i  =  3) 

*  superscript  used  to  indicate  the  value  of  any  variable  evaluated  when  a  point  (a|,  Oq)  in  the 
diagram  of  temperature  gradients  is  on  RS 

**  superscript  used  to  indicate  the  value  of  any  variable  evaluated  when  a  point  (aj,  Oq)  in  the 
diagram  of  temperature  gradients  is  on  R}* 
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Growth  Condition  of  an  Ice  Layer  in 
Freezing  Soils  Under  Applied  Loads 
2.  Analysis 


YOSHISUKE  NAKANO  AND  KAZUO  TAKEDA 


INTRODUCTION 

In  this  report  we  will  consider  the  one-directional  steady  growth  of  an  ice  layer.  Let  die  freezing 
process  advance  from  the  top  down  and  the  coordinate  x  be  positive  upwards,  with  its  origin  fixed 
at  some  point  in  the  unfrozen  part  of  the  soil.  A  freezing  soil  in  this  problem  may  be  considered  to 
consist  of  three  parts:  the  unfrozen  part  the  frozen  fringe  Rj  and  the  ice  layer  R2,  as  shown  in  Figure 
1.  The  physical  properties  of  parts  jRq  and  are  well  understood  but  our  knowledge  on  the  physical 
properties  and  the  dynamic  behavior  of  part  Rj  does  not  appear  sudicient  for  engineering  applica¬ 
tions. 

It  has  been  shown  empirically  (Radd  and  Oertle  1973,  Takashi  et  al.  1981)  that  there  is  a  unique 
temperature  T*  *  at  ni  for  a  given  pressvire  of  ice  P21  at  nj  and  a  given  pressure  of  water  Pj  in  Ro  when 
an  existing  ice  layer  neither  grows  nor  melts  and  the  mass  flux  of  water /i  in  R]  vanishes.  This  temper¬ 
ature  Ti*  at  the  phase  equilibrium  of  water  is  given  as 

a=P2i-Pi  =  -YTr,  if  /i  =  0  (1) 


where  y  is  a  constant  with  the  value  of  1.12  MPa  ®C-’,  and 
a  and  P21  are  often  referred  to  as  the  effective  pressiue 
and  the  overburden  pressvire,  respectively.  Equaticm  1  is 
often  called  the  generalized  Clausius-Clapeyrcm  equa- 
tiim,  which  is  attributed  to  Edlefsen  and  Anderson  (1943). 

Konrad  and  Morgenstem  (1980, 1981, 1982)  empiri¬ 
cally  found  tiiat  the  rate  of  water  intake at  the  forma¬ 
tion  of  the  final  ice  lens  is  proportional  to  the  average 
temperature  gradient  (T\  in  die  frozen  fringe.  This  may 
be  written  as 


/io  =  -SPo(r).  (2a) 

where  a  prime  denotes  differentiaticm  with  respect  to  x. 
The  positive  proportionality  factor  SPg  is  termed  the  seg¬ 
regation  potentiaL  which  is  a  property  of  a  given  soil. 
Konrad  and  Morgenstem  also  fovind  empirically  ttiat 
SPg  is  a  decreasing  functicm  of  bod)  die  applied  pressure 
o  and  the  suction  of  water  (-Pjg)  at  ng.  We  will  write  this 
dependence  as 


Figure  1.  Schematic  of  a  steadily  growing 
ice  layer  in  a  freezing  soil. 


erally  depend  on  the  tempeiahire  T  and  the  composititm  of  the  soil.  We  will  describe  sudi  functional 
dependence  of  (i  =:  1^)  as 


fC,  =  K2(T,p) 

(4d) 

K2  =  K2(T,p) 

(4e) 

where  p  symbolically  denotes  the  composition  of  the  soil  that  is  uniquely  determined  by  the  bulk 
densities  of  unfrozen  water  Pi,  ice  P2  and  soil  minerals  P3. 

Nakano  (1990)  obtained  the  exact  mathematical  solution  to  the  problem  of  a  steadily  growing  ice 
layer  based  on  the  model  M].  Analyzing  the  behavior  of  this  solution,  we  showed  that  Mj  is 
consistent  with  eq  1  (Nakano  1990)  and  eq  3b  (Nakano  and  Takeda  1991).  We  also  showed  (Nakano 
and  Takeda  1991)  that  Mi  can  accurately  describe  the  steady  growth  condition  of  an  ke  layer  under 
negligible  applied  pressures. 

The  steady  growth  condition  of  an  ice  layer  with 
or  without  applied  pressure  is  the  tegitm  Rebound* 
edbyacurveRs  and  a  straight  line  Rs*  in  the  dia¬ 
gram  of  temperature  gradients  as  shown  in  Figure 
3.  The  region  R,  is  defined  as 

(lc,/l!o)a,  >  Oo  >  ki  (fcb  +  (5) 

where  ki  and  kg  are  the  thermal  conductivities  of  R2 
and  Rq,  respectively,  Oq  is  the  absolute  value  of  the 
temperature  gradient  at  Mq  and  O]  is  the  limiting 
value  of  the  temperature  gradient  at  ni  in  R2  de¬ 
fined  as 

oi  =  -  lim  ^(a^.  (6) 

X-MIl 

xmX2 

In  eq  5,  L  is  the  latent  heat  of  fusion  of  water,  b  is  a 
function  of  the  thickness  5  of  Ri  defined  by  eq  73c  and  73e  in  Nakano  (1990),  and  R21  the  limiting 
value  of  K2  X  approaches  ni  when  x  is  in  Ri  and  an  asterisk  denotes  that  K21  is  the  value  of  K21  when 
a  point  (ai,Oo)  is  on  Rs  in  the  diagram  of  temperature  gradients. 

In  Figure  3  we  will  refer  to  the  region  as  R,„  where  Oq  >  (ki  /  ko)ai  ,/io  <  0  and  an  ice  layer  is  melting. 
The  botmdary  Rs  *  is  given  as 

ao  =  (ki  /ko )  tti  on  RJ  * .  (7a) 

An  existing  ice  layer  neither  grows  no  melts, /iq  varushes  and  eq  1  holds  true  on  RJ  * .  The  botmdary 
Rs  is  given  as 

<Xo  =  ki(ko  +  o*'  •  (7b) 

It  is  easy  to  see  from  eq  5  that  the  steady  growth  condition  of  a  given  soil  is  uiuquely  determined  by 
O]  and  Oq.  Nakano  (1990)  showed  that  all  physical  variables  such  as  f^g,  Tj,  5,  etc.,  are  also  uniquely 
determined  by  O]  and  Og  for  given  hydraulic  conditions,  and  applied  pressure  o.  The  hydraulic  con¬ 
dition  in  our  tests  is  specified  by  the  distance  Sg  between  ng  and  n  where  the  pressure  Pi„  of  water 
is  kept  at  the  atmospheric  pressure.  Therefore,  any  point  in  R,  is  urtiquely  specified  by  otj  and  Og  for 
given  and  o.  We  will  write  this  as 


Figure  3.  Temperature  gradients  Oj  and  Og. 
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R,  =  /?,(o„ao;  6o,Pi„,a). 


(8a) 


Since  Oj  and  Oo  are  related  by  eq  7b  on  R? ,  any  point  on  R2  is  uniquely  specified  by  either  a,  or  Og 
for  givoi  Sg,  Pin  and  o  as 

RJ*RJ(oq;  oaPi,-4  (8b) 

The  main  objective  of  this  work  is  to  show  that  M)  is  ccxisistent  with  the  data  under  various  applied 
pressures  that  were  presented  in  Part  I.  We  will  also  show  that  Mj  is  consistent  with  the  reported 
empirical  equations  such  as  eq  2a-c  and  3b-c.  Furfoermoie,  we  will  evaluate  the  concept  of  seg¬ 
regation  potential  introduced  by  Konrad  and  Morgenstem  (1980, 1981, 1982)  based  on  Mj  and 
experimental  data. 

PROPERTIES  OF  Ml 

Treating  a  given  soil  as  a  mixture  of  water  in  liquid  phase  B},  ice  Bi  and  soil  minerals  Bj,  Nakano 
(1990)  obtained  the  exact  mathematical  solution  to  the  problem  of  a  steadily  growing  ice  layer  based 
on  foe  model  Mi  under  foe  following  assumptions.  The  density  of  each  constituent  remains  constant, 
foe  dry  density  of  Rg  renuins  constant,  foe  part  Rg  is  kept  saturated  with  water  at  all  times  and  foe 
pressure  of  water  Pi^  at  some  boundary  n  fixed  in  Rg  remains  constant.  The  thermed  conductivity  k(x) 
in  R|  is  assumed  to  be  a  nondecreasing  linear  function  of  x  given  as 


lc(x)  =  *gIl+Tl(X-ng)], 

»!  >  X  2  Hg 

(9a) 

Tl  =  (*bi-*o)/(S*fl) 

(9b) 

Um  k-ko\^  ki 

xm  Ki 

(9c) 

8  =  «,  -  «g. 

(9d) 

The  temperature  T  in  Ri  satisfies  foe  equation  given  as  (Nakano,  1990) 


k{x)T'  -  Cifio  T  =  -  kgOg . 


(10a) 


The  solution  of  eq  10a  is  approximately  given  as 


T(ad  =  -otoJ(x  -  no)  +  i  (po  -  q)  (x  -  no)*j 

(10b) 

Ti=-0tg<l(8) 

(10c) 

r(nt)  =  -agiK8) 

(lOd) 

where  ^(n  f )  is  foe  limiting  value  of  ^(x)  as  x  approadies  while  x  is  in  Ri,  and  a,  b  and  Pg  are  defined 
as 


(11a) 


«(5)  =  8  +  (l/2)(pg-ii)52+... 
iK5)  =  (l  +  il8)-»lH-po5+...] 
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(11b) 


Po® 


(11c) 


We  derived  the  following  equation  of  heat  balance  in  jR|  given  as 

kjOi  *  IhoOo  +  (1-  “  C2Ti)/io  (12) 

The  mass  flux  of  water  satisfies  the  equations  (Nakano  1990)  given  as 

Pio  =  l’in~l(/io/l^  +  Pol^  (13a) 

r"l  rn\ 

P2i  =  Pio-/io  Xr‘dr-j  Xr’Xjrdx  (13b) 

Jno  Jhq 

where  P,o  =  P,(no),  Pi„  =  Pi(n) 
n  s  some  point  in  Pg 
Xg  =  the  hydraulic  conductivity  in  Pg 
Pa  =  gravity  term,  density  di  x  gravitational  acceleration 


8o  =  iig-n^0.  (13c) 

In  order  to  reduce  eq  13b  to  a  simpler  form,  we  introduced  (Nakano  and  Takeda  1991)  the  following 
two  dimensionless  quantities: 

♦o(T)  =  r'[  (Xio/Xi)(X2/X2o)dT  (14a) 

Jo 

♦i(T)  =  r»[  (Kxo/Ki){k/ko)dT  (14b) 

Jo 

where  X^g  and  Kjo  are  the  limiting  values  of  Xj  and  X2,  respectively,  as  x  approaches  ng  while  x  is  in 
P|.  We  obtained  (Nakano  and  Takeda  1991)  die  following  equations  given  as 

Xg  =  Xjg  (15a) 

Y=X2o/Xg  (15b) 

♦02  =  1,  if/2g  =  0  (15c) 

-Ti  =  (<T  +  5oXg*/io )  /  f  (16a) 

1=  Y«oi-X^^j4n  (16b) 

0=  P21-Pl,  (16c) 

where  ^oi/  and  y  are  defined  as 

♦w  =  ^b(r,)  (17a) 
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♦ll  -  ♦l(A) 

(17b) 

y=/io/«o- 

(17c) 

Equation  15a  holds  true  because  the  composition  of  the  freezing  soil  is  continuous  at  %  and  eq  15b 
and  c  follow  from  eq  4b. 

For  the  sake  of  convenience  we  will  reduce  eq  16a  to  a  form  similar  to  eq  3b.  First,  we  will  introduce 
die  following  four  dimensionless  quantities; 


[O,  77*  =  0 

|(-  r  f  [l  -  K2(yKi)1  dT,  7t*  <  0 


(18a) 


[(KoK2)(KiK2o)i 


dZ 


ri  =  7t* 


Ti<7i* 


(18b) 


0)2  = 


[0,  77*  =  0 

1 

(-77*)  [l-kmoKif\dZ  77*  <0 

i  Jo 

Ih  Ti  =  77* 


fTi 

-  AT"*  kKo(koKi)  dT,  Ti  <  Ti* 

i 


(18c) 


(18d) 


Using  eq  18a-d,  we  will  write  and  as 


Ti^oi  =  '(7*(1  +  ©o)  -  ©1  AT  (19a) 

=  77*(1  +  ©2)  -  ©3  AT.  (19b) 

According  to  Mj  the  mass  flux  of  water  /ig  in  a  neighborhood  of  nj  in  J?]  is  given  as  (Nakano  and 
Takedal991) 

fio  =  -K„P;(fti>-)-^bK2iao  (20) 

where  Ku  and  Pi(nf  )are  the  limiting  values  of  Kj  and  Pi,  respectively,  as  x  approaches  ni  while  x 
is  in  P}.  We  will  rewrite  eq  20  as 

y  =  ^f2i<Pb  (21a) 


<|H,=  6-K„P,'(n,+)(oofC2ir’ 


(21b) 
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We  found  (Nakano  and  Takeda  1991)  that  is  p<»itive  in  R,  and  vanishes  on  Ks*;  namely 
Pi'(n?-)>0  in  R,  (22a) 

P,'(ni*‘)  =  0  on  Rj.  (22b) 

From  eq  22a,  22b  and  21b  we  obtain 

q)o=  b  on  Rs  (23a) 

b  xpo  >  0  in  Rs  (23b) 

<po  =  0  on  R?*.  (23c) 

Using  eq  21a,  we  will  reduce  eq  16b  to 

I  =  lK^oi-ri%M  (24a) 

where  Yt  is  defined  as 

Yi=K2iKm.  (24b) 

Substituting  ^  and  in  eq  24a  by  eq  19a  and  b  and  using  eq  1,  we  will  reduce  eq  24a  to 

-Ti/  =  «t[1  +  Wo  -  Yi<po(l  +  W2)]  +  Y[ttti- YKPotaslAT.  (25a) 

Combining  eq  25a  with  16a,  we  obtain 

5o^o‘Vio  +  [Yi<Pb(l  +  ®2)  -<»o]  o  =  y[wi  -  Yi  900)3]  AT.  (25b) 

Now  we  will  reduce  eq  25b  to  a  form  similar  to  eq  3b  as 

AT  =  So  +  fli  fw  (25c) 

where  ao  and  a\  are  defined  as 

2o  =  [yi  90(1  +  wj)  -  Wo]  0(19 1  )■'  (25d) 

«l  =  «o(Y^^o9l^'  (25e) 

9i  =  0)i -Yi9oO)3-  (25f) 

We  will  examine  eq  25c  for  a  special  case  where  fig  is  very  small.  It  follows  from  eq  21a  fiuit  <po 
vanishes  as  fig  vanishes.  We  find  from  eq  18a  that  tOg  vanishes  because  of  eq  4b  as  /lo  vanishes.  It 
follows  from  eq  1  that  T^  approaches  Ti* ;  hence,  o)}  approaches  one  as  fig  approaches  zero.  HeiK«, 
ao  approaches  zero  and  a\  approaches  8b(Y^S))~*  os  fig  approaches  zero.  Therefore,  when  fig  is  very 
snu^  eq  25c  may  be  approximately  given  as 

AT^SiAo.  (26) 
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Since  /iq  is  the  growth  rate  of  the  ice  layer,  eq  26  states  that  the  growttt  rate  of  the  ice  layer  is 
proportional  to  the  degree  of  supercooling  and  that  the  rate  coefficient  (ai)"*  depends  on  the 
hydraulic  properties  of  Rq  and  Rj,  namely  the  availability  of  imfrozen  water.  Equation  26  is  ccmsistent 
with  the  theory  of  crystal  growth  in  supercooled  liquid  (Chalmers  1964). 

We  will  detoe  the  thicknesses  6,  of  Rn  and  8**  of  Rjo  (Fig.  2)  as 

8e  =  «i-«io  (27a) 

8**  =  «,o-«o.  (27b) 

Then  the  thickness  5  of  Ri  is  obviously  given  as 


5  =  5“+ 8*. 

Using  eq  10c,  we  obtain 


(27c) 


=-aofl(8‘*). 


(27d) 


Using  eq  10c  and  11a,  we  obtain 


ATsOode 


(27e) 


where  a,  is  defined  as 

a,  =  fl(8)-fl(8**) 

(27f) 

=  8e  +  l  (Po-tl)(8e  +  28**)8e+.... 

(27g) 

Using  eq  27e,  we  will  reduce  eq  26  to 

«e  =  3iy 

(28a) 

where  a\  may  be  written  as 

2i  =  SoiyKoT^  (oai  -  Km  a}3y)r\  (28b) 

When/io  is  small,  y  and  8,  are  also  small  and  eq  28a  and  b  ate  reduced  to 


8e=Siy 

(29a) 

2i  =  8o(YKoa)i)'^. 

(29b) 

From  eq  29a  and  b  we  tod  that  the  thickness  of  fite  essential  frozen  fringe  is  proportional  to  y  and 
dtat  the  essential  frozen  fringe  vanishes  as/^g  vanishes.  In  odier  words,  when/jg  is  very  small,  we  may 
state  that  the  appearance  of  an  essential  frozen  fringe  is  induced  by  die  flow  of  unfrozen  water 
regardless  of  o.  This  implies  that  an  essential  frozen  fringe  appears  only  under  a  dynamic  condition. 

MODEL  Ml  AND  SEGREGATION  POTENTIAL 

We  will  show  below  that  Mj  is  consistent  with  eq  2a,  which  was  found  empirically  by  Kcmrad  and 
Morgenstem  (1980, 1981)  and  was  confirmed  by  Ishizaki  and  Nishio  (1985).  In  a  typical  experiment 
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Figure  4.  Schematic  0/  trajectories  oi(t)  that 
approximately  descrilK  the  condition  of  freez¬ 
ing  at  the  formation  of  the  final  ice  lens. 

0 

by  Konrad  and  Morgenstem  (1981),  the  temperature  field  in  the  s)^tem  changes  rapidly  at  the  start 
of  the  experiment.  However,  as  time  elapses,  the  rate  of  the  change  slows  down  so  that  the  transient 
freezing  may  be  accurately  approximated  by  a  series  of  successive  steady  states.  Hence,  the  later  part 
of  the  experiment  can  be  approximately  represented  by  trajectory  1  in  Figure  4,  consisting  of  points 
a(f)  =  {ai(0,  ao(t)l  for  ti^t^  to- 

A  point  a(fo)  is  in  where  frozen  soil  without  any  visible  ice  layer  grows.  As  decreases  and 
Ofl  increases,  the  trajectory  approaches  the  vicinity  of  a  point  a(f  j)  in  R^.  As  we  described  previously 
(Takeda  and  Nakano 1990)  the  pattern  of  ice-rich  frozen  soil  grown  in  this  vicinity  evidently  depend^ 
on  the  soil  type  and  the  magnitude  of  (or  Oq).  The  restilts  of  tests  on  Kanto  loam,  for  instance,  clearly 
indicate  that  the  pattern  of  rhythmic  ice  banctog  is  formed  at  the  small  values  of  Ui  while  soil  particles 
or  small  aggregates  of  soil  particles  are  evenly  dispersed  at  the  greater  values  of  (X]. 

When  a(f)  reaches  the  point  a{ti)  on  RS,  the  final  ice  layer  emerges.  While  a{t)  moves  toward  the 
point  0(12)  on  Rs*,  the  growth  of  the  final  ice  layer  continues  with  the  decreasing  growth  rate  until 
(x(f)  reaches  the  point  ol{1^  on  Rs  *  where  the  ice  layer  stops  growing.  It  should  be  noted  that  a  line 
of  constant /u  is  nearly  parallel  to  RS  *  because  of  eq  12.  From  eq  20  and  22b  we  obtain  on  RS 

fio  =  y*ao=  Kjibao  on  RS.  (30) 

It  follows  from  eq30  that  the  water  intake  flux,  fjo,  at  the  formation  of  the  final  ice  layer  is  proportioruil 
to  the  temperature  gradient,  to(y  at  .  Comparing  eq  30  with  2a,  we  find  drat  SPg  and  (T),  in  eq  2a 
correspond  to  Kji  ~  ^  respectively.  Since  the  temperature  gradient  in  Ri  does  not  vary 

significantly,  the  segregation  potential  SPg  is  nothing  but  K21  (die  limiting  value  of  the  transport 
function  K2asx  approaches  while  x  is  in  Rx  at  die  formation  of  the  final  ice  layer),  when  a  point 
a(t)  is  on  Rs*  in  dw  diagram  of  temperature  gradients,  namely 

SPo=Kii^yb-\  (31) 

We  have  shown  diat  Mx  is  consistent  with  eq  2a.  It  is  clear  from  eq  20  that  eq  30  holds  true  on  RS  but 

does  not  hold  in  R,  because  of  eq  22a.  In  other  words,  the  value  of  y  defined  by  eq  17c  is  equal  to  bKn 

on  RS .  However,  Ae  value  of  y  in  R,  depends  on  a  specific  trajectory  a(t).  For  instance,  on  trajectory 
1  in  Figure  4  die  value  of  y  decreases  from  bK^i  as  0((t)  moves  toward  the  point  odt^  from  the  point 
a(f  x)  and  vanishes  at  die  point  adjtfj.  On  diis  trajectory  fig  decreases  widi  the  increasing  Og.  However, 
it  is  easy  to  see  diat/xg  decreases  with  the  decreasing  0^  on  trajectory  2  for  ^2  >  f  >  fi-  Therefore,  Mx 
is  also  consistent  wifo  die  empirical  finding  by  Ishizaki  and  Nishio  (1985)  that  the  vai  ue  of  y.  varies 
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widely  and  that  /jg  may  either  increase  or  decrease  with  the  increasing  Og  depending  on  a  given 
specific  trajectory.  This  finding  by  Ishizaki  and  Nishio  (1985)  was  empirically  confirmed  (Nakano 
and  Takeda  1991)  when  o  =  0. 

According  to  our  definition  K21  is  the  value  of  K21  evaluated  on  Rs*.  Since  any  point  on  R/  is 
uniquely  specified  by  eq  8b,  we  may  write  X21  ^ 

K2*  =  K2*i(aQ;  So,  P^,,  a).  (32) 

We  studied  (Nakano  and  Takeda  1991)  the  dependence  of  X21  on  otg  for  a  special  case  where  Sg  =  2.0 
cm.  Pin =0.1  MPa  and  o = 0  for  three  types  of  soils.  It  was  found  that  K21  is  nearly  constant  for  a  p 
soU  if  cxg  is  greater  than  2.0  (°C  cm-i).  However,  K21  tends  to  increase  with  decreasing  Og  in  the . 
of  Og  with  2.0  (“C  cm-*)  >  Og  >  0  for  the  two  types  of  soils,  Tomakomai  silt  and  Fujinomori 
Unfortunately,  we  were  imable  to  confirm  behavior  similar  to  this  of  X21  for  Kanto  loam  because  of 
a  lack  of  data.  Using  the  additional  data  newly  obtained,  we  will  show  such  behavior  for  Kanto  loam 
below. 

Konrad  and  Morgenstem  (1980, 1981, 1982)  empirically  found  eq  2b  that  is  equivalent  to  the 
following  equation  given  as 

K2*i  =  K2*,(pro,  a)  (33) 

where  an  asterisk  for  Pig  is  used  to  emphasize  that  the  value  of  Pig  is  evaluated  when  a  point  a(l)  is 
on  R  * .  Since  their  hydraulic  conditions  were  not  specified  in  the  same  marmer  as  in  oiir  experiments, 
we  will  reduce  eq  33  to  the  form  appropriate  to  our  system.  In  our  system  Pig  is  given  by  eq  13a. 
Hence,  we  obtain 

Pio  =  Pi«-[(^o/Ko)  +  po]8o  (34) 

where /ig  isthevalueof/igOnRg*  and  is  uniquely  determined  by  ag(eq  8b)  if  5g,  Pi  „  and  o  are  given. 
Therefore,  Pig  in  eq  33  can  be  replaced  by  Og,  80  and  Pi„  so  that  eq  33  is  reduced  to  eq  32.  We  have 
shown  that  M]  is  consistent  with  eq  2b,  which  was  found  empirically  by  Konrad  and  Morgenstem 
(1980,1981,1982). 

Konrad  and  Morgenstem  (1981)  empirically 
fotmd  that  SPg  is  a  monototucally  decreasing 
function  of  -Pig  when  a  =  0.  Since  Pig  repre¬ 
sents  the  combined  effects  of  Og,  and  Pi^  in 
order  to  find  the  effect  of  the  temperature  gra¬ 
dient,  we  plotted  the  data  of  SPg  vs.  (-T% 
obtained  by  Konrad  and  Morgenstem  (1981)  in  SPg 

Figure  5,  where  the  number  assigned  to  each 
data  point  corresponds  to  the  test  number  of 
their  E  series  experiments.  The  data  points  E4 
through  7  were  obtained  for  a  single  layer  of 
Devon  silt  imder  varioris  temperature  condi- 
timis.  These  data  points  clearly  indicate  that  SPg 
tends  to  increase  with  the  decreasing  tempera¬ 
ture  gradient.  This  tendency  is  consistent  with 
our  empirical  findings.  Tests  E8  and  E9  were  (-T') , 

both  two-layer  systems  in  which  the  hydraulic  figure  5.  Values  o/SPg  [g(cm  ^  dh*/  »s.  the  values 
conductivitiesKooftheunfrozenbottomlayers  of  the  average  temperature  gradient  (-T ),(•<:  cm-i) 

were,  respectively,  higher  and  lower  than  that  obtained  by  Konrad  and  Morgenstem  (1981). 
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of  the  vinfrozen  Devon  silt.  Since  die  temperature  gradients  in  diese  two  tests  were  not  equal,  the 
effects  of  hydraulic  conditions  alone  on  SPq  are  difficult  to  assess. 

RESULTS  OF  DATA  ANALYSIS 
Steady  growth  condition 

We  will  excunine  the  validity  of  the  model  Mj  under  die  applied  pressure  by  using  the  experimen¬ 
tal  data  presented  in  Part  I  (Takeda  and  Nakano  1993).  We  found  empirically  that  the  steady  growth 
region  Rjla)  under  a  given  applied  pressure  a  is  approximately  described  as 

au  =  Aat,  kiko^  >  A>  S(a)  (35) 

where  Ou  =  absolute  value  of  the  temperature  gradient  near  no  in 

Of  =  absolute  value  of  the  temperature  gradient  near  n^  in  R2 
S  =  property  of  a  given  soil  that  depends  on  the  applied  pressure  a. 

The  temperature  profiles  in  the  frozen  and  the  unfrozen  parts  are  not  exacdy  linear  when  die 
steady  growth  of  an  ice  layer  is  taking  place  because  of  the  convective  heat  transport.  However,  the 
amount  of  heat  transported  by  convection  is  much  less  than  that  transported  by  conduction.  The 
difference  between  ol^  (or  oCf)  and  Oq  (or  a^)  is  negligibly  small  as  shown  empirically  and  dieoretically 
in  Nakano  and  Takeda  (1991).  Therefore,  eq  35  is  nearly  equivalent  to 

ao  =  Aai,  -kiko  >  A>  S(a).  (36) 

We  will  no  longer  discriminate  a„  (or  o^)  from  Oo  (or  Uj)  in  the  following  discussion. 

According  to  Mi  the  steady  growth  region  R,(o)  under  a  given  applied  pressure  a  is  given  by  eq 
5 .  Using  y*,  we  will  reduce  eq  5  to  a  form  similar  to  eq  36  as 

ao  =  Aai,  kik^  >  A>  ki(fco  +  Ly*)  .  (37) 

Comparing  eq  37  with  36,  we  find  that  Mi  is  consistent  with  the  experimental  data  if  die  following 
relation  holds: 

S(o)  Sr  jfci[jto  +  Ly*(<T)r-  (38a) 

We  will  define  y  *  as 

S(o)=jti[jto  +  Ly*(o)r.  (38b) 

Then,  it  is  easy  to  see  that  eq  38a  is  equivalent  to  the  following  relation: 

y*(o)^y*(o).  (38c) 

We  will  examine  the  validity  of  eq  38c  below. 

Foragivenothe  value  of  y*(o)canbe  calculated  by  usingdiecalciilated  value  of/iobased  on  eidier 
die  measured  rate  of  heave  r  or  die  measured  rate  of  water  intake,  and  die  measured  value  of  Oq  at 
each  data  point  (O],  0^  on  R^.  The  calculated  values  of  y*(o)  are  plotted  vs.  Oq  with  a  beit^  a 
parameter  in  Figure  6.  The  mass  flux  of  water  fu  decreases  by  die  order  of  1(H  as  o  increases  from 
zero  to  195  kPa.  As  die  result  die  accuracy  of  measuring/10  tends  to  decrease  widi  increasing  o  and 
die  variability  of  data  points  becomes  mote  pronotmced  with  increasing  o,  as  shown  in  Hgure  6. 


Figure  6.  Values  of  y*  [g(cm  “C  dh’i  os.  oto  under  various  applied 
pressures  a  (kPa). 


Because  of  such  variability  and  the  limited  numbers  of  data  under  a = 390 kPa,  we  decided  not  to  use 
the  data  taken  imder  a  =  390  kPa  in  our  analysis.  From  Hgure  6  we  find  the  general  trend  that  y* 
decreases  with  increasing  a  and  that  y*  increases  with  the  decreasing  Oq  in  the  range  of  Oq  <  2.0  "C 
cm~>.  The  latter  trend  was  also  observed  (Nakano  and  Takeda  1991)  in  the  experiments  with 
Tomakomai  silt  and  Fujinomari  clay  tmder  null  applied  pressure. 

We  calculated  the  values  of  y*(o)  from  the  values  of  S(a)  that  were  presented  in  Table  2  of  Part 
I  (Takeda  and  Nakano  1993).  We  also  calculated  the  average  of  y*(o)  over  all  the  data  points  obtained 
for  eacho.  The  values  ofy*(o)  and  die  average  values  ya*(o)ofy*(<T)  for  eachoare  presented  in  Table 
1.  It  is  clear  from  Table  1  that  eq  38c  does  not  hold  for  every  o,  particularly  for  greater  values  of  a. 
However,  die  average  values  ya*(o)  do  not  difier  significantly  from  dioseof  y*(o).  We  may  conclude 
that  the  model  is  consistent  with  the  experimental  data  regardless  of  oand  that  die  steady  growth 
region  of  an  ice  layer  under  various  applied  pressures  can  be  described  by  eq  37. 

In  order  to  find  the  dependence  of  ya*on  a,  we  plotted  ya*in  the  logarithmic  scale  against  a  in 
Figure  7.  It  is  clear  from  ^e  figure  that  ya*is  a  decreasing  function  of  o.  Assuming  diat  b  is  nearly 
equal  to  one,  we  may  conclude  that  fC2i  (or  SPq)  is  a  decreasing  function  of  a,  which  was  found 
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Figure  7.  Average  values  yj  Igfcm  d)-*J  os,  o 
(kPalufhereasoUdlineis  tluempirically  determined 
relationship  for  Devon  silt  obtained  by  Konrad  and 
Morgenstem  (1982). 


Table  1.  Calculated  valucsy*  (glcm^CdHl 
and  the  average  measured  values  y,  under 
various  applied  pressures  o  (kPa). 


AygUtd  gratmts  (a) 


0.0 

8.12 

16.2 

48.7 

97.8 

195 

r 

2J4 

2.12 

1.14 

1.12 

0.72 

0.56 

y* 

2.86 

1.98 

1.61 

0.97 

0.45 

035 
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empirically  by  Kcmrad  and  Morgenstem  (1982).  Their  data  with  Devon  silt  were  presented  in  the 
functional  form  given  as 


SPo  =  a2exp(-fl30).  (39) 

When  we  use  the  same  tmitsof  SPq  and  a  as  those  ofya*and  o  in  Figure  7,  the  values  of  the  constants 
a2  and  a3  are  1.04  and  8.95  x  10^,  respectively,  and  eq  39  is  presented  by  the  straight  line  in  Figure 
7.  Because  of  the  limited  number  of  data  points  it  is  not  certain  that  our  data  can  be  presented  in  the 
same  functional  form  as  eq  39.  However,  the  important  point  is  that  K21  is  a  decreasing  function  of 
a.  The  reason  for  such  dependence  will  be  discussed  below. 

Dependence  of  y*  on  7? 

Combining  eq  30  with  eq  32,  we  obtain 

y •  =  bKa  (aa  8a  Pi„,  o)  •  (40a) 

For  a  special  case  such  as  our  experiments  where  Sq  and  Pi„  are  specified,  we  may  reduce  eq  40a  to: 

y*  =  8X2*1  (aa  o).  (40b) 

On  the  other  hand  X2*  is  the  value  of  X21  when  a  point  (ttj,  o^,)  is  on  P** .  From  eq  4e  we  obtain 

X2*  =  ffaVTi*,  p)  (40c) 

where  T*  is  the  temperature  at  n,  when  a  point  (a,,  c^)  is  on  R^ .  It  is  clear  from  eq  40b  and  c  that  T* 
and  the  composition  p  generally  depend  on  Oq  and  o.  We  will  study  empirically  the  relationship 
between  y*  and  T*  below. 

Using  the  set  of  data  obtained  imder  the  applied  pressure  of  48.7  kPa  as  an  example,  we  will 
describe  how  we  obtained  the  empiricedly  determined  value  of  T*  from  the  data.  The  results  of  our 
data  aiuilysis  are  presented  in  Table  2,  where  «i  is  the  observed  location  of  the  interface  between  R, 
and  Ry  while  ng  is  the  location  of  the  0°C  isotherm  calculated  by  using  the  measiured  temperature 
profile  in  Rq.  The  values  of  8  in  the  table  are  calculated  simply  from  eq  9d  and  vary  between  0.91  and 
1.5  trun.  We  have  found  that  the  value  of  8  increases  with  foe  increasing  a  and  foat  foe  maximum 
value  of  8  in  foe  range  of  o  ^  195  kPa  is  4.8  mm  imder  foe  condition  of  Og  =  0.80  and  a  =  195  kPa. 

The  value  of  T*  can  be  calculated  from  either  foe  measured  temperature  profile  in  Rg  or  foat  in  R2. 
As  we  discussed  in  Part  I  (Takeda  and  Nakano  1993),  foe  temperature  measurements  in  Rg  are  more 
accurate  than  foose  in  R2.  Therefore,  it  is  desirable  to  determine  T*  from  foe  profile  in  Rg.  However, 
the  thermal  conductivity  k(x)inRi  is  unknown  because  foecompositionofRi  is  unknown.  According 
to  foe  model  M^,  is  given  by  eq  10c  and  11a.  Hence,  when  foe  variation  of  Ic  in  R^  is  small,  T*  is 
nearly  equal  to  T,o  defined  as 


Table  2.  Summary  of  data  analysis  with  o  s  48.7  kPa. 


Oo 

rCcnr>) 

fio 

(g  crn^iH) 

y* 

[g<cm<C 

"1 

<cm) 

"o 

(cm) 

6 

(cm> 

“Tio 

rc) 

-T„ 

rc) 

A* 

-Ti 

rci 

1 

0.642 

0.968 

151 

0.67 

055 

0.12 

0.076 

0.174 

0.125 

2 

1.29 

1.59 

123 

0.25 

0.099 

0.15 

0.192 

ai98 

0.195 

3 

1.62 

1.60 

0.988 

028 

0.13 

0.15 

0239 

0211 

0.225 

4 

2.76 

2.13 

0.770 

0.18 

0.086 

0.094 

0.260 

0.245 

0.255 

5 

329 

265 

0.604 

0.20 

0.080 

0.12 

0578 

0.320 

0549 

6 

5.85 

2.97 

0507 

0.13 

0.039 

0.091 

0535 

0.467 

0.501 
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(41a) 


It  is  clear  from  eq  10c  and  11a  that  T*  is  accurately  approximated  also  by  Tm  when  6  is  very  small. 
The  calculated  v^ues  of  are  listed  in  Table  2.  When  o  is  negligible,  Nakano  and  Takeda  (1991) 
found  that  k{x)  tends  to  increase  with  x  in  Ri-  Therefore,  Tjq  would  be  a  lower  boimd  of  T*,  namely 

Tf  2  Tio.  (41b) 

We  also  calculated  the  value  of  T*  from  the  measured  profile  in  Ri-  The  calculated  values,  which  are 
referred  to  as  Tn,  are  listed  in  Table  2.  We  find  from  Table  2  that  Tjo  tends  to  be  less  than  Tjj.  A  ten¬ 
dency  similar  to  this  was  also  found  in  all  other  cases  of  different  applied  pressures.  Under  these 
circumstances  we  decided  to  choose  the  average  of  T ,q  and  Tj  j  to  be  the  empirically  determined  value 
fr  of  Tf  defined  as 

fi*  =0.5(Tio  +  Tn).  (42) 


The  values  of  Ti  are  listed  in  Table  2.  ^ 

The  values  of  y*  are  plotted  against  -Ti  with  the  logarithmic  scale  under  various  applied  pres¬ 
sures  a  in  Figure  8.  Despite  some  scatter,  it  is  clear  that  the  relation  between  y*  and  Ti  is  nearly  one 
to  one.  The  solid  line  in  Figure  8  is  the  best  linear  approximation  to  the  data  points  given  as 

y*=K*=/^“  (43a) 

lK2d(fo/T)^  r<fo  (43b) 


wheie  b  is  assumed  to  be  one,  Kjo  is  the  limiting  value  of  as  x  approaches  Mq  while  x  is  in  Rj  and 
is  equal  to  1.98  x  10*g(cm  d  "C)-!  (Nakano  and  Takeda  1991),  To  =  - 1.5  x  10-*°C  and  =  1.039.  As 
we  showed  (Nakano  and  Takeda  1991),  K20  satisfies  the  equation  given  as 


where  Kq  is  the  hydraulic  conductivity  in  Rg.  It  is  easy  to  see  that  y*  becomes  infmite  as  T  approaches 
zero  in  eq  43b.  Although  eq  43b  is  the  best  approximation  to  the  data  points,  eq  43a  is  needed  to  fit 
the  data  points  in  a  neighborhood  of  T  =  0°C. 


Figured.  Valuesofy*[g(cm  ‘CdHJvs.  the  temper¬ 
ature  -Ti  (“C)  under  various  applied  pressures  a 
(kPa). 
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According  to  eq  40c,  y*  generally  depends  on  Tf  and  the  composition  p.  However?  we  have  found 
empirically  that  y*  depends  mainly  on  T* .  This  implies  that  the  composition  in  a  neighborhood  of 

in  I?]  is  not  significantly  affected  by  a  and  Oq  (or  fig).  Since  T i  varies  tetween  0  and  -1 .0°C  in  Figure 
8,  we  may  conclude  l  .at  the  composition  of  the  essential  frozen  fringe  Rn  depends  mainly  on  the 
temperature  regardless  of  a  and/jQ.  There  is  another  interpretation  of  Figure  8:  that  the  transport 
function  K2  does  not  depend  on  the  composition.  Recently  Nakano  and  Tice  (1990)  found  empirically 
that  K2  in  unsaturated  frozen  clay  strongly  depends  on  the  composition,  particularly  the  content  of 
ice.  Their  empirical  finding  supports  the  former  interpretation. 

Figure  8  shows  that  the  range  of  T\  for  a  given  o  shifts  toward  the  lower  temperature  as  a  increas¬ 
es.  The  segregation  potential  K21  evidently  depends  primarily  on  the  temperature  T*  at  ttj  and  is  an 
increasing  function  of  T*  because  the  applied  pressure  in  the  range  of  a  ^  195  kPa  does  not  affect  sig¬ 
nificantly  the  composition  of  the  essential  frozen  fringe  Ru.  This  is  the  reason  why  the  segregation 
potential  K21  is  generally  a  decreasing  function  of  o. 

Dependence  of  T*  on  fio 

The  values  of-Ti  are  plotted  agai^t  f^g  under  various  applied  pressures  o  in  Figure  9.  Figure  9 
shows  that  the  relationship  between  T 1  and  fig  is  approximately  linear  for  a  given  a,  which  is  con¬ 
sistent  with  the  empirical  relation  (eq  3b  and  c)  found  by  Ishizaki  and  Nishio  (1985).  It  is  clear  horn 
Figure  9  that  the  constants  ag  and  aj  strongly  depend  on  a.  An  important  question  is  whether  we  can 
describe  the  behavior  of  data  points  in  Figure  9  by  eq  25c  derived  based  on  Mj.  We  are  not  able  to 
show  that  eq  25c  is  consistent  with  the  data  because  we  have  no  data  on  the  transport  function  K,. 
However,  we  will  show  below  that  eq  25c  is  consistent  with  the  data  if  the  function  Ki  is  properly 
chosen. 

The  model  Mj  is  defined  as  the  frozen  fringe  where  ice  may  exist  but  does  not  grow,  and  the  mass 
flux  of  water /j  is  given  by  eq  4a  with  the  condition  of  eq  4b  and  c.  When  o  >  0,  the  essential  frozen 
fringe  Rjj  vanishes  as/i  vanishes  but  Rjo  does  not.  From  eq  4b  we  obtain: 


X2/Ri  =  y  inR,oif/,  =  0.  (45) 

When  the  steady  growth  of  an  ice  layer  occurs,/,  remains  constant  at /,o  throughout  R^,  R^  and  R,,. 
An  important  question  arises:  whether  or  not  eq  45  holds  true  when /jq  does  not  vanish.  In  other 


0.8 


0.6 


0.4 

Figures.  Values  of CO  vs.  the  mass 
fluxcfwaterfiglg  (cnr^  rf-*)/  under  vari¬ 
ous  applied  pressures  a  (kPa). 

02 


0  2  4  6  8 

fTo 


15 


words,  is  dw  composition  of  Rjo  significantly  affected  by /^o?  As  we  described  above,  our  data  indi¬ 
cate  that  neither  a  nor /,o  significantly  aff^  the  composition  of  the  essential  frozoi  fringe.  There¬ 
fore,  it  is  probable  that  the  efiect  of on  the  composition  of  Rjo  ^  negligible.  We  will  assunne  foat 
eq  45  holds  true  regardless  of /jg,  namely 

Xj/ Xj  =  y  in  Rig.  (46) 

When  eq  46  holds  true,  the  dimensionless  quantity  defined  by  eq  18a  vanishes.  Hence,  eq  25d  is 

reduced  to 


So^ynpoU  +a)2)c(Y<ih)  *• 


(47) 


Now  we  will  calculate  T*  as  a  function  of  by  eq  25c  as  follows.  Since  T*  is  the  value  of  Ti  on  R/, 
tpig  is  equal  to  5  by  eq  23a.  We  will  assume  that  6  =  1,  or  equivalently  k{x)  -kgUiRi.  We  also  assume 
that  X2(T)  is  equal  to  X2i(T)  given  by  eq  43a  and  b.  The  value  of  T*  *  is  calculated  by  eq  1  for  a  given 
o.  Horiguchi  and  Miller  (1983)  found  empirically  that  the  transport  functions  X^fT)  of  various  soils 
can  accurately  be  represented  in  the  same  functiorud  form  as  eq  43b.  We  will  assume  that  Xi(T)  is 
given  as 


X,(7)  =  /^" 

IxoCfo/rf 


To<  T^O 
To 


(48) 


where  Kg  is  Ibe  hydraulic  conductivity  in  Rg  and  is  1.77  x  1(P  g(cm  d  MPa)-i  (Nakano  and  Takeda 
1991).  The  value  of  To  is  the  same  as  used  in  eq  43a  and  b.  A  coirstant  bj  is  an  unknown  parameter 
to  be  determined. 

As  we  described  in  Part  I  (Takeda  and  Nakano  1993),  the  applied  pressure  o  affects  the  void  ratio 
e  of  a  specimen.  Although  the  variation  of  e  itself  is  negligibly  small,  the  hydraulic  conductivity  Kg 
may  be  affected  significantly.  Therefore,  we  determined  empirically  foe  relationship  between  Xg  and 
a  given  as 


Xg(o)  =  1.77  X  KP  0-0“*'  (49) 

where  foe  imits  of  Xg  and  a  are  g(cm  d  MPa)-*  and  kPa,  respectively.  The  value  of  Xg  is  reduced  to 
about  one-half  according  to  eq  49  when  o  is  increased  from  zero  to  195  kPa.  The  functioruil  form  of 
eq  49  is  consistent  with  foe  data  obtained  by  Fukushima  and  Ishii  (1986).  In  our  calculations  of  T* 
for  a  >  0  we  used  Xg(o)  given  by  eq  49  instead  of  foe  value  of  Xg  at  o  =  0. 

Now  we  can  calculate  T* (fto)  wifo  6,  being  a  parameter.  Calculating  T,(^,g)  in  foe  wide  range  of 
6],  we  find  that  foe  calculated  curves  T*  (fio)  fit  foe  data  well  if  6]  is  about  one-half  of  62.  The  calculated 
curves  wifo  61 = 0.52  are  presented  in  Figure  10  together  wifo  foe  data.  If  6]  is  decreas^  (or  increased) 
firom  this  value,  foen  foe  gradients  of  foese  curves,  d(-  Ti)/df^o,  increase  (or  decrease).  We  have 
shown  that  eq  25c  is  consistent  wifo  foe  data  if  foe  function  X]  is  given  by  eq  48  wifo  6}  =  0.52. 

DISCUSSION  AND  CONCLUSIONS 

Many  models  of  frost  heave  have  been  proposed  in  foe  past  (Nakano  1990).  However,  foe  model 
proposed  by  Konrad  and  Morgenstem  (1982)  is  rate  of  few  that  were  built  on  an  empirical  base.  Their 
segregation  potential  theory  was  easily  adapted  to  solve  engineering  problems  in  foe  past.  As  our 
quantitative  understanding  on  foe  subject  is  increased,  their  model  can  be  improved  or  refined 
without  sacrificing  its  easy  adaptability  to  engineering  problems. 

Konrad  and  Morgenstem  (1981)  proposed  an  equation  similar  to  eq  4a  where  the  mass  flux  of 
water/i  is  given  as  foe  sum  of  two  terms,  namely,  a  pressure-related  term  and  a  temperature  related 
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FigurelO.Vttluesof-Ti  (’C)vs.thefluxf*Q[g 
(cmr^d-^)]tt^e  solid  curves  are  predicted  re¬ 
lationships  betroeen  these  two  variables  by  eq 
25c  under  various  applied  pressures  a  (kPa). 


term.  However,  curiously  enough,  dtey  dropped  the  pressure  term  in  their  later  publications.  Since 
the  pressure  term  in  eq  4a  is  generally  negative  according  to  M],  the  omission  of  diis  term  leads  to 
overestimating/i.  Therefore,  their  model  certainly  predicts  an  upper  bound  of  frost  heave  as  they 
claim  (Konrad  and  Morgenstem 1982).  However,  some  serious  criticisms  against  their  model  cannot 
be  refuted  unless  die  pressure  term  is  restored,  as  we  will  explain  below. 

When  the  pressure  term  is  neglected,  it  is  dear  from  eq  4a  diat/i  is  nonnegative.  This  is  die  reason 
why  Konrad  and  Morgenstem  (1982)  could  not  provide  a  satisfactory  explanation  for  the  expulsion 
of  water  from  freezing  soils.  Talcashi  et  al.  (1978)  conducted  a  series  of  frost  heave  tests  in  which  the 
temperature  in  die  unfrozen  part  Kq  was  kept  omstant  at  02-03‘’C  higher  than  the  freezing  point  of 
specimens.  In  other  words,  die  positive  temperature  term  of  eq  4a  was  kept  small  in  their  tests. 

The  absolute  value  of  the  negative  pressure  term  of  eq  4a  is  small  when  the  applied  pressure  o  is 
small.  Henoe,/i  can  be  positive  when  o  is  small.  However,  the  pressure  term  decreases  with  the  in- 
creasingoand/iVanishesatcertainvalue8ofobecausethetwotermsofeq4acancelout.AsaincreaS' 
es  beyond  this  value,/]  becomes  negative;  that  is,  the  expulsion  of  water  from  freezing  soils  takes 
place.  Takashi  et  al.  (1978)  found  empirically  ndut  we  described  above. 

Anodier  serious  criticism  of  die  segregation  potential  theory  was  raised  by  Ishizaki  and  Nishio 
(1985)  diat  the  value  of  y,  defined  by  eq  3a  is  constant  stricdy  at  die  instant  when  the  final  ioe  lens 
emerges,  but  exc^t  for  this  instant,  y,  varies  widely  during  die  growdi  period  of  the  final  ice  lens. 
As  we  explained  earlier,  die  pressure  term  of  eq  4a  vanishes  at  the  formaticm  of  die  final  ice  lens,  but 
the  negative  pressure  term  varies  depending  on  a  specific  trajectory  in  the  diagram  of  temperature 
gradients  when  the  final  ice  lens  is  growing.  The  empirical  finding  by  Ishizaki  and  hfishio  (1985)  can 
explained  if  die  pressure  term  is  restored.  Assuming  that  die  temperature  Ti  (o)  at  die  formation 
of  the  final  ice  lens  depends  mainly  on  the  property  of  a  given  soil  alone,  Konrad  and  Morgenstem 
(1982)  termed  T* (o)  as  the  "segregation  freezing  temperature."  However,  we  have  found  empirical¬ 
ly  (Fig.  9)  that  it  (<1)  depends  strongly  cm  the  mass  flux  ho¬ 
using  the  data  obtained  in  Part  I  (Takeda  and  Nakano  1993),  we  evaluated  the  accuracy  of  Mj .  We 
found  that  die  predicted  steady  growth  ccmdition  of  an  ice  layer  under  various  applied  pressures  is 
in  gcxxl  agreement  with  diat  found  empirically.  We  also  found  diat  M]  is  consistent  with  die  data 
obtained  by  Konrad  and  Morgenstem  (1980,  IS^l,  1982)  that  were  used  to  support  their  segregatum 
potential  theory.  Their  mediod  of  frost  heave  predicticm  is  a  sound  and  useful  tcx>l  for  engineering 
problems  diat  consistendy  provides  an  upper  bound  of  frost  heave.  However,  the  accuracy  of  their 
method  can  be  significantly  improved  and  scene  of  the  serious  criti-cians  against  it  can  be  refuted 
if  die  pressure  term  in  die  equaticm  of  water  flow  is  restored  as  we  discussed  above. 
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